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ABS TRACT . 0 ft 

1 

The later stages of stellar evolution, when the internal 

temperature is above 10 9 °K, are affected by neutrino processes. 

Under these conditions the neutrino emission occurs at a rate 

higher than the photon emission rate. It is found that the 
homology relation (Density) ~ (Temperature)^ breaks down completely 
and that in general the density rises much faster than the third 
power of temperature. In regions where nuclear burning commences 
<C 12 - 0 16 or 0 16 - 0 16 reactions) the interior stellar structure 
is determined by the condition: (Nuclear Energy Production Rate) 

+ (Neutrino Energy Dissipation Rate) ■ 0. Such curves have nega- 
tive power dependence on temperature (p ~ T~ 12 ) or are isothermal 
depending on the density. The rapid removal of energy from the 
core of the star results in a short time scale for these stages 
and the star continues to contract. The stars evolve (independently 
of their mass) to a very high central density (~ 10 9 "" 10 g/cm 3 ) at 
a rather low central temperature (~ 3-4 x 10 9 °K) . It may be 
concluded that the star does not reach the full iron-helium phase cha 
as previously suggested (Fowler and Hoyle, 1964) . Stellar collapse 
will be induced at this density and temperature either by p-transi- 
tion or by a small amount of iron- helium transition. 
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INTRODUCTION 

The problem of supernova collapse has acquired considerable 

interest in recent years, in connection with general relativity 
and element synthesis. In most literature conclusions are based 

on extremely simplified structures . for presupernova stars. The 
most often quoted presupernova model is a polytropic model of 
polytropic index n = 3, and the subsequent contraction to a dynamical 
collapse is assumed tp be homologous. Xn the later stages of 
stellar evolution when the temperature exceeds 10 K, it has been 
shown that neutrino processes will dissipate stellar energy rapidly, 
with a time scale for evolution of the order of a few thousand 
years or less, this time scale is shorter than the time scale for 
redistribution of radiation energy in the star. On the other band 
in the homologous contraction hypothesis it is implicitly assumed 
that the time scale for energy redistribution is much shorter than 
the time scale for stellar evolution. Therefore, homologous 
contraction hypothesis cannot be applied to stellar evolution 
problems when T ~ 10 9 °K. It is the purpose of this paper to 

discuss the structure of presupernova stars without assuming the 
validity of the homologous contraction hypothesis. 

THE BASIC EQUATIONS 
The stellar structure equations are 



(1) 

(2) 

(3) 

(4) 


The notations used are the standard notations (e.g. those used 

by Schwarzschild (1958)), 

We now discuss the influence of neutrino processes on stellar 
structure. The cross section for interaction of neutrinos with 
matter is around 10 ^ cm 2 , in contrast to that for photons which 
is around 10" 20 -* 1(T 25 cm 2 . The mean free path for neutrinos 
in lead is around one light year! Only under extreme density 
conditions (p £ 10 11 g/cm 3 ) must one consider the absorption of 
neutrinos by stellar matter (Chiu 1964) . Neutrino opacities in 
dense matter have been considered by Bahcall (1964) and Bahcall 
and Frautschi (1964). Since the star does not absorb neutrinos, 
the only place they can enter the stellar structure equations is 
through equation (3) which describes the energy balance of the 
star. In equation (3) S> is the energy production term, and we have 


(5) 


e * e gr + e n + (du/dt) , , 

where 6g r is the rate of gravitational energy release (in ergs/g-sec) , 

o n ulta »«()• or nunion^ anevyy tJoletiaa, and (dU/dt) is the rate 

\ . 

energy is dissipated by neutrino emission. 

we shall now estimate the rate of photon energy transfer and 

the rate of neutrino emission inside a star. The photon opacity 

*ph stellar matter is always greater than n e = 0.19 cm /g 

(the opacity for electron scattering) • The mean free path of photons 

is roughly (1 /h s ) in g/cm^, and- hence is always shorter than 5 g/cm^. 

The luminosity of a star apart from a numerical factor which is 

not very different from unity is roughly 

t (Photon Energy Content of the Star) 

(r/ c) • (R/X) ' (6) 

where R is the radius of the star and c is the velocity of light, 
and X. is the photon mean free path (in cm) . The relaxation time 

for a star to be cooled by photon emission Tp^ is therefore 

. v _ (Photon' Energy Content of the Star) c= ( 7 ) 

T ph " (L) -■ cX * , 

Using a simple model in which the temperature T and the radius R 
are related by the following condition (which is derivable from the 
virial theorem) 

- (Gravitational Energy) «= (GM°/R ) = (3/2)R TM 

9 

' • ** (Thermal Energy), i 


we find 


' sec 


( 9 ) 


t “ 2xlO l i t 
ph 0 




The corresponding relaxation time via pair annihilation neutrino 

emission t is 
V 

_ . _ = 1.5 x 10 4 exp[ll .9 /Tq]Tq (g/M) 2 , (Tg « 6).(10) 

T v - (dU/dt) y y 

At T 9 = 1, t v = 2xl0 7 sec and Tph = 2X10 11 sec. Hence we can 
conclude that at T 9 ~ 1 the photon energy transport process cannot 
be important in the center of stars. 

Soon after the temperature exceeds ~ 10^ °K the rate of energy 
dissipation by neutrinos will dominate over the photon processes. 

Thus there exist only narrow temperature-density regions where 

one must consider both photon and neutrino processes simultaneously. 

Outside this region one can either completely neglect the photon 

energy transfer process (except for luminosity calculation) or 

completely neglect the neutrino processes. We shall define a 

star to be a photon star if the major part of the energy loss of 

the star is through surface emission, and a neutrino star if the 

major part of the energy loss is through neutrino emission. 

The idea of a neutrino star is only an abstraction: In real 

3 5 o 

stars the surface temperature is in between 10 -10 K, thus xn 
the outer region of a star photon processes must dominate. However, 
the neutrino active core will evolve rapidly before the outer 
region can effectively affect the core's structure (e.g. by the 
addition of mass to the core, or the addition or removal of energy 


from the core) . Thus, the concept of a neutrino star will be helpful 
in obtaining the evolutionary properties of a star in its later 
stages. 

Since we neglect L as compared with - (dU/dt) dM r (where 
the integration extends over the whole star) we can also assert 
that in the neutrino active region (dL/dM) « -(dU/dt) , In the 
first approximation to equation (3) we replace (dL/dM) by zero. 
Equation (3) then becomes 

e - - T(dS/dt) M + & n + (dU/dt) «= 0, (11) 

where S is the entropy per unit mass.*' Partial differentiation of 
S with respect to t is taken at a shell of constant mass My. 
(Convection cannot occur in neutrino stars, because the temperature 
gradient never exceeds the adiabatic temperature gradient.) 

Effectively equation (4) may be ignored in calculating the 

. 

mechanical structure of the neutrino star in regions where neutrino 
emission is strong. That is to say, the structure of the star is 
not effected by the inclusion of equation (4) . Once the structure 
is obtained one can calculate L from equation (4) . The equations for 
the stellar structure of a neutrino star are given by equations (1) 
and (2) and equation (11). It can be shown that if p(r) and T(r) 
are given at t = t 0 , the structure of the star at a later time is 
uniquely determined. The proof is quite elementary, though tedious. 

\ 


*It can easily be demonstrated that £ gr = - T(dS/dt) . 
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EVOLUTION OF A NEUTRINO STAR 
WITHOUT NUCLEAR ENERGY GENERATION 

Tho energy sources of stars are either gravitational or nuclear 
in origin. Nuclear reactions occur at well defined temperatures. 


and when nuclear energy is not available the star must contract 
to supply the energy lost via neutrino emission or photon emission. 

*7 A 

At lower temperatures (T £ 10 K) it was found that stars 
contract in a simple fashion. Under the assumption that the temper 
ature and density relations of a star are such that one can write 


p = K (U(l/n)) (12) 
whore K is a constant depending on the mass of the star and n is 
the polytropic index, one can simplify equations (1) and (2) by 
the following substitutions 


. r 0 5 > 

where r Q and \ are also constants. 
Lane-Emden equation 


(13) 

(14) 

The resulting equation is the 


l _d_ 




(15) 


However, a family of solutions is derivalbe from one existing 
solution for equation (15) ; If 9(§) is a solution, then it can 
be shown that 

A (2/(n-l) ) Q(A§) . ( 16 ) 

is also a solution (this is known as the homology hypothesis) . 
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Further, when one applies this hypothesis to evolving stars, one 
finds that inside a star, at a given shell mass M^, the relation 
of the density p and the temperature T in the seu ree of 

P * T 3 . ■ (17) 

Hence, if the polytropic index n is a constant in the course of 
stellar evolution, then equation (17) is valid and the structure 
of a star can be characterized by the function 8(§) ; the evolution 
will introduce only a change in the scale. (For proof, see for 
example, Eddington (1930) or Fowler and Hoyle (1964).) 

This simple result works remarkably well for many types of 
stars. One implicit assumption involved in the application of 
the homology hypothesis to stellar evolution (this is called the 
homologous contraction hypothesis) is the following: Energy can 

br redistributed inside stars at a rate much faster than it can 
be produced by gravitational contraction. 

Fowler and Hoyle (1960, 1964) have applied the homologous 
contaction hypothesis to stars in their late evolutionary stages. 

Under certain assumptions they concluded that a star of mass around 
30 Mq will collapse when the center reaches a temperature of ~ 6x10^ °K 
and a density of 2x10^ g/cm^. A separate estimate by Chiu (1964) 
and by Chiu and Fuller (1962) using a similar hypothesis showed 
that long before this temperature and density is reached a collapse 
will take place via neutrino processes. 



-10- ' 

We now examine whether the homologous contraction hypothesis 
can be applied to neutrino stars. In the absence of £ n equation (11) 
becomes 

- («JU/OU) << 0) ‘ 


• ® “ (c /p)T^ for nondegenerate stars 
v 0 

where To ~ 3, e =4.3x 10 15 , and K = 9. For a nondegenerate 
9 v o 

gas, S is given by the following equation 

S «= Rglnt^p" 1 ) + Constant , (19) 

where a is a slowly varying constant; a = 1.5 for nonrelativistic 
matter (Tg « 1) and a = 3 for relativistic matter (Tg 6 ) ; at 

Tq b 2, ft W 4 • / • 

We now assume that at each shell mass p and T evolve according 

to the following relation 

(p/p 0 ) = (T/T )B- (20) 

u o> 

where a , T , and ft are constants. The homologous contraction 
K o o 

hypothesis will require that 0 = 3 throughout the star. Substi- 
tuting equations (19) and (20) into equation (18) one obtains 


(dT 8 /dt) = Y T 0 K “6 


( 21 ) 


where 


S 0 T o 


Pq R g (8 ” a ^ 


Equation (21) can be solved exactly to give 


t 8 = (t 8 °)/(i - (t/t 0 )C V(K-a-i) ] 


where 



x is the relaxation time for the neutrino energy loss (equation 
v 

(10) ) . Thus for a homologously contracting neutrino star, a 
collapse occurs at a time t = t Q . If we use 3 = 3, K » 9, a = 2.7, 
then t Q = (1/45 ) t v . The free fall time for a 10 M Q star is reached 
T 9 ~ 3 (Chiu 1964) . 

The collapse can be delayed to higher temperatures if the 

value of 3 is close to K - 1 = 8 . One exact calculation on the 

evolution of neutrino stars was recently done by Chiu and Salpeter 

(1965) . This calculation confirms our guess that the value of 3 

is close to 8 for the central region of the star. The stellar 

mass is chosen to be 10 M 0 and the initial configuration is chosen 

3/2 , 

to be that of a polytrope of n = 1.5(p « T ) . No approximation 
was made to the equation of state which included contributions 
from electrons and positrons, radiation, and heavy nuclei (iron). 

In the log p - log T plane the structure of the star is represented 
initially be a straight line of slope 1.5. If the star evolves 
homologously, then this line will be displaced without distortion 
along a line of slope 3 (Figure 1) . The actual evolution of the 
star is very complicated, as indicated in Figure 2. Since the 
neutrino energy loss rate at high densities is from the plasma 
neutrino process, the rate increases with increasing density, and 
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a positive temperature gradient develops throughout the inner 
1/3 of the mass of the star. 

Xn ■kvhin Vn hnv«i dAmoncHsir erfcrac} -hbr*^ A.* nufilatir naSitAi&BH 

are neglected then the homologous contraction hypothesis cannot 
bo applied to neutrino starr. , for control conditions evolve according 
to p ~ T n where 6 £ n £ 9 . 

EVOLUTION OF NEUTRINO STARS WITH NUCLEAR REACTIONS 
The time scale t v for a star with no energy source other than 
the gravitational energy to remain at a temperature T is given by 

• « E/(-clU/dT) . (25) 

In reality the' time scale is shortened by a factor of around 4. 

On the other hand nuclear energy (from 0^ - 0-^ reaction and the 
C^ 2 - 0^ reaction) is fairly abundant. The time scale T n for 
nuclear burning at a temperature such that 

e n + (dU/dt) <=0 (26) 

is several times (~ 10) greater than t hut still much less 
than the photon diffusion time. Because nuclear reactions have 
much steeper temperature dependence (~ T^®) than that of the 
neutrino rate (~ T^) a negligible amount of nuclear fuel will 
be consumed below the temperature defined by equation (26) . 

Once nuclear burning starts, in that part of the star , 
gravitational contraction can take place only along curves given 
by the condition 


/ 
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e gr + l n + < dU / dt) “ 0 


(27) 


Since, in 6g r , p and T appear in logarithms (equation (19)), the 
eonaitlen 

/dlj\ •%# 

* & n ■ ‘ (Constant ) -1 , ( 28 ) 


determines the temperature and density relation according to which 
the star can evolve. 

For the C 12 - 0 16 and - O 3 ' 6 nuclear reactions, Fowler 

and Hoyle- (1964) have given the following rates:* 

lo g & 12-I6 “ 49 * 7 + lo %° pX 16 X 12 “ 12 / 3 ) 10 ^ T 0 

- (46.30/T 0 V3) (1 + o.080T e )Va i (29) 
lo 9 e l6-16 15 55 * 7 + lo Sio P X 16 X 16 " (2.3)log . 0 T 9 


' - (59.04/T 8 ^ 3 ) (1 + 0.080T o )^ 3 , (30) 

0 3 

where x is the concentration of element a. At p — 10 g/cm we 
a 

find that the temperature defined by equation (26) is roughly 

16 13 

T 9 = 2.45 for the 0 16 - 0 16 reaction and Tg = 1.25 for the 0 - C 

reaction r at these two temperatures one can approximate £-^6— 16 an< 3 
^12-16 ky simpler formulae 


£ 16-16 = 8 xi 0 ” 6 P T o 3 ° ergs/g-sec , 
C 12-16 = i .32xl0 3 pT a 3 1,8 ergs/g-sec 


(31) 

(32) 


Equation (28) thus gives 


*At present calculations by Salpeter and Deinzer (1964) , Stothers 
(1964) and Divine (1964) seem to indicate that C 12 - C 12 and C 12 -0 16 
reactions do not take plkce, since carbon will not be produced by 
helium burning in massive stars. 



-14- 


t 

OxlO“ 6 pT 0 3a + (4.3xiO‘Vp)T 0 ° - (Constant) *T 

(0 16 - o 16 Burning) , (33) 

' 4 

• 1 .32xlG 3 pT 9 a 1,8 + (4.8xl0 i0 /p)T o a cxp[- 2(T o /T 0 )] 

= (Constant) -7 (C 12 “ ° 16 Burning). (34) 

In the log p - log T plane equations (33) and (34) are lines of 

negative slope, which we shall call the n-v curves. The particular 

curves with the constant replaced by zero will be called the n-v-zero 

curves. We have plotted the n-v-zero curves for 0^ - 0 ^ burning 
12 16 

and C - 0 burning (Figure 3) . 

To what degree can the structure of a neutrino star in the 

nuclear burning stage itself deviate from the n-v-zero curves? 

The nuclear energy rate has a very sharp temperature dependence. 

A deviation in temperature T by a factor of 4 % from the curve 
will , cause nuclear energy production rates to vary by more than ... 
a factor of 3 while the neutrino rate only varies from by about 
40 % i Imagine that the star is on the lower temperature side; 
then a gravitational contraction will raise the temperature of 
the star until nuclear energy production takes over the gravitational 
energy production and in this part of the star contraction will 
start. If the star is on the higher temperature side, then the 
excess energy produced by nuclear reactions will expand the star 
until the temperature is low enough so that an energy balance is 
reached. Thus the star is always pushed toward the n-v-zero line. 
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Because of the extremely rapid cooling by neutrino emission, 
there is a positive temperature gradient in the core of the star 

and the maximum temperature, and eo nuclear burning,, ocoure in a 
shell surrounding the core. Thus the core continues to collapse 

with p ~ T® until it reaches the n-v-zero line, and nuclear burning 
begins. The nuclear burning shell evolves along the n-v-zero line. 
And, because the time scale of evolution is much shorter than the 

photon diffusion time, the energy from the nuclear burning is not 

transferred to the envelope, which also continues to contract. 

Thus the structure of the star approaches the n-v-zero line. When 
the central density approaches 4 x 10^ g/cm 3 for 0^ - 0^ burning 
or 4 x 10 6 g/cm^ for C 3-2 - O 16 burning, the plasma neutrino process 
becomes important. The plasma neutrino loss rate Q t can be approxi- 
mated to an accuracy of better than 50 % at x £ 1 by 

“ 6.04 + 3 log T 9 + log p 6 , ; (35) 

, where if pg » 1, x is given by 

x 0*237(1 + 0.6 4 130 6 2/3 r< P 6 *T 9 ; ; (30) 

The n-v-zero curve becomes a vertical line and the temperature is 
independent of the density and is a constant. For the - O^-O 
burning process we have 

■ l°ff - O.I 38 - (l/30) log(x^x b ) , Tq ~ I. 38 , (37) 

Time, scale ~ 10^ sec , (38) 

Thus during nuclear burning the core evolves according to the 
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condition & n + (aU/dt) = 0 and a dense core is formed independent 
of the mass of the star . Only an exact calculation will tell 

us what the final density will be« Until 1 in 

completely consumed the temperature of the star cannot increase. 

The time scales for evolution via neutrino emission in a 
realistic star are indicated on Figure 3. These differ drastically 
from those predicted on the basis of simple polytropes. A neutrino 
star, it seems, evolves quickly (in 10* 7 sec) at a temperature of 
around 1.2 x 10 9 °K (C 12 - 0 16 burning) or 2.5 x 10 9 °X (0 16 - 0 16 
burning) , its density increasing rapidly. The time scale then 

becomes longer (~ 10® - 10® sec for 0^ - 0^ burning, and ~ 10^ sec 
for C - 0 burning) . It is not known if the density will reach 
a value of ~ 10 19 g/cm 2 . 

At densities around 10 9 g/cm 3 @-transitions and pycnonuclear 
reactions (Cameron 1959) will cause elements to approach an 
equilibrium state quickly. Tsuruta (1964) has considered the 
equilibrium of elements in this density-temperatrue regime. This 
approach to nuclear equilibrium or, at higher density the forma- 
tion of neutrons, or even a small amount of Fe disintegrates to He, 
may be the trigger of the supernova collapse. 

On the basis of evolutionary properties of photon stars Burbidge, 
Burbidge, Fowler, and Hoyle (1957) , and later in a more elaborate 
version Fowler and Hoyle (1960, 1964) , have argued that the 


transition from iron to helium at a density of around IQ* 7 g/cm® 

and a temperature of 7 x 10 K is sufficient to cause a star of 

mass around 30 Mg to collapse. However, it seems that this density 
and temperature cannot be reached simultaneously in neutrino stars; 

the structure of neutrino stars is so drastically different from 
that of photon stars that their argument is invalidated. It is 
possible that the final collapse state will be at a very high 
density (~ 10®*"^ g/cm^) and at a much lower temperature (~ 3-4x10® °K) 

CONCLUSION 

We have found that the classical homology hypothesis p 
is completely inapplicable to the evolution of stars whose tempera- 
tures are over 10 K, where neutrino processes dominate over 
photon processes - such stars we have called neutrino stars. 

Instead, neutrino stars during gravitational contraction evolve 
with p c ~ T ®, while during nuclear burning they evolve along the 
n-v-zero curves which lead any star (independent of its mass) to 
a region of very high density. After the exhaustion of nuclear 
fuel the central temperature will increase again. The final collapse 

of all neutrino stars into supernova will probably occur at a density 

4 — i o •a 9 o 

of around 10 g/cm and at a temperature of around 3-4 x 10 K. 

The mechanism for collapse is thought to be the approach to equi- 
librium of elements (at high density) via 0 -transitions or pycno- 


nuclear reactions 
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We have not included the electron conduction opacity. 

• 

At p ~ 10 9 g/cm 3 the opacity due to electron conduction is ex- 
ceedingly small, the positive temperature gradient may thus be 
replaced by an isothermal gas sphere . 
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LEGENDS 

i 

Figure 1. The evolution of homologously contracting polytrope 

model stars with a polytropic index n “ 1*5 » All parts ef the star 
evolve according to the law p « T^. For a 10 Mg star the center 

will reach the Fe^ - He^ region at a density p ~ 10^ and T — 7x10^ K. 
While the concept of homologous contractions works reasonably well 
for photon stars,, it breaks down completely for neutrino stars as 
in Figure 2 . 

Figure 2. The evolution of neutrino stars. The structures at 
various times are shown* The highest density point on each curve ( 
is the center of the star. Although t_ is given to 7 decimal 
places only the time differences between stages are useful numbers. 

The evolution is very different from that for a photon star (Figure 1) . 
In this model the collapse condition is obtained at p ~ 10 9 g/cm 3 
while the star is still composed of iron. This model demonstrates 
the complete break down of the concept of homologous contraction 
in neutrino stars . 

Figure 3. n-v— zero curves for 0^ — and 0^ - 0^ reactions, 
x is the concentration of the element _a and a curve with x^ *= 0.2 
is also shown for comparison. Along the n-v-zero curves the nuclear 
energy release rate is equal to the neutrino rate. The bending 
at higher densities occurs because the plasma neutrino rate becomes 
important and the plasma neutrino rate has a different density ~ 
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temperature dependence from that of annihilation neutrinos. Num- 
bers on the graphs refer to (log <r n ** log x fl ) (r n is the time scale 

in seconds for nuclear burning at that point) . The 60 % « He 

curve is also shown. Dotted line is the evolution 

track of a model star of 30 M Q studied by Fowler and Hoyle (1964) . 
Actual stars must evolve along the n-v-zero lines. The hyphenated 
line (- — - — - — - — )' separates regions in which plasma 
or annihilation neutrino processes dominate as indicated. 


) 






(g/cm 5 ) 








